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Abstract 

We study the dynamics of first-order phase transition in the early Universe 
when it was 10 — 50/is old with quarks and gluons condensing into hadrons. 
We look at how the Universe evolved through the phase transition in small 
as well as large super cooling scenario, specifically exploring the formation 
of quark nuggets and their possible survival. The nucleation of the hadron 
phase introduces new distance scales in the Universe, which we estimate 
along with the hadron fraction, temperature, nucleation time etc. ft is 
of interest to explore whether there is a relic signature of this transition 
in the form of quark nuggets which might be identified with the recently 
observed dark objects in our galactic halo and account for the Dark Matter 
in the Universe at present. 
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It is well known that a phase transition from quark gluon plasma to confined 
hadronic matter must have occurred at some point in the evolution of the early 
Universe, typically at around 10 — 50/is after the Big Bang. This leads to 
an exciting possibility of the formation of quark nuggets through the cosmic 
separation of phases H . As the temperature of the Universe falls below the the 
critical temperature T c of the phase transition, the quark gluon plasma super 
cools and the transition proceeds through the bubble nuclcation of the hadron 
phase. The typical distance between the nucleated bubbles introduces a new 
distance scale to the Universe which depends critically on the super cooling 
that takes place. This bubble nucleation of the hadron phase is governed by the 
formation of critical bubbles, that is bubbles bigger than the critical size start 
to expand and coalesce till the entire Universe is filled with the hadron phase. 
As the hadronic bubbles expand, they heat the surrounding plasma, shutting 
off further nucleation and the two phases coexist in thermal equilibrium. The 
hadron phase expands driving the deconfined quark phase into small regions of 
space and it may happen that the process stops after the quarks reach sufficiently 
high density to provide enough pressure to balance the surface tension and 
the pressure of the hadron phase. The quark matter trapped in these regions 
constitute the quark nuggets. The number of particles trapped in the quark 
nugget, its size and formation time are dependent sensitively on the degree 
of super cooling. The duration of the phase transition also depends on the 
expansion of the Universe and on other parameters like the bag pressure B and 
the surface tension a. 

The quark nuggets formed in the small super cooling scenario are in a hot 
environment around the critical temperature T c and are susceptible to evapo- 
ration from the surface Q and to boiling through subsequent hadronic bubble 
nucleation inside the nuggets [^| . However in the large super cooling scenario we 
have the interesting possibility of these nuggets forming at a much lower temper- 
ature than T c due to the long time of phase transition and consequent expansion 
of the Universe. Alcock and Farhi Q have shown that the quark nuggets with 
baryon numbers < 10 52 and mass < 1O~ 5 M0 are unlikely to survive evapora- 
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tion of hadrons from the surface. Boiling was shown to be even more efficient 
mechanism of nugget destruction by Alcock and Olinto ^ . If these nuggets are 
formed at around lOOMeV, they cannot have more baryons than are allowed in 
the horizon by the standard model of Cosmology, i.e. ~ 10 49 ( 10Q J IeV )~ 3 and 
would not survive till the present epoch. We will see that for small super cooling 
this is the likely fate of the nuggets. These results were somewhat modified by 
Madsen et.al. GJ by taking into account the flavor equilibrium near the nugget 
surface for the case of evaporation and by considering the effect of interactions 
in the hadronic gas in a relativistic mean field model described by Walecka 
for the case of boiling || . They suggested that the quark nuggets with large 
baryon number allowed by the causality limit may after all be able to survive 
from the early Universe. In the large super cooling scenario the time at which 
these nuggets are formed can be quite late along with a much lower formation 
temperature due to the expansion of the Universe. Such nuggets can have a 
baryon number content of > 10 52 and being at a much lower temperature than 
T c (< 0.1T C ) will easily survive till the present epoch. The number of baryons 
in the horizon (of size ~ 2t) can also be large due to long time it takes for the 
bubbles to meet and form the nuggets. 

There have been recent observations by gravitational micro lensing jf| of 
dark objects in our galactic halo having masses of about 0.01 — 1 solar mass. 
If these objects have to be identified with quark nuggets, they could only have 
been formed at a time later than the time when the Universe cooled through T c , 
i.e. later than 50-100 /is after the Big Bang. Such a possibility of the nuggets 
forming at a temperature ~ O.lMeV, implying a high degree of super cooling 
and strongly first order phase transition was recently investigated by Cotting- 
ham et. al. jjj in the Lee- Wick model [8J. Their investigation showed that the 
time of formation and the baryon content of these nuggets are essentially deter- 
mined by the rate at which the hadron bubbles nucleate. However there was still 
the question of reheating due to the expansion of the bubbles of hadron phase 
which raises the temperature towards T c . These studies have also been carried 
out by taking interactions into account in both the phases and by incorporating 
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the effects of curvature || energy in the calculations. 

Studies of quark-hadron phase transitions in the early Universe |Tc| | , in heavy 
ion collision JTl] ] and in high density nuclear matter have been done previously by 
looking in detail at the dynamics of the phase transition. Kapusta et. al. have 
applied a recently computed nucleation rate |IlJ to a first order phase transition 
in a set of rate equations to study the time evolution of a quark-gluon plasma as 
it converts to hadronic matter in heavy ion collisions. Based on Bjorken hydro- 
dynamics and on current parameter values, they find the transition generates 
30 percent extra entropy and also a time delay of ~ 11/m/c in completion of 
the transition. Kajantie and Kurki-Suonio [jl0| studied how the early Universe 
during the quark-hadron phase transition evolved through the mixed phase in a 
scenario with small initial super cooling. Fuller et.al.|T(| have also studied the 
dynamics of the Universe during the constant temperature coexistence epoch. 

In this letter we calculate the nucleation rate and quantitatively study what 
happens when the temperature drops to T c . This nucleation rate is used to solve 
a set of rate equations to study the time evolution of the quark-gluon phase as 
it converts to hadronic matter in an expanding Universe. A novel feature of 
this method is that reheating of the plasma during phase transition is included 
in the calculations and all relevant quantities can be evaluated as a function of 
time. 

Bubble Nucleation:- When the early Universe as a quark-gluon thermody- 
namic system cools through the critical temperature T c , energetically, the new 
phase remains unfavorable as there is free energy associated with surface of sep- 
aration between the phases. Small volumes of new phase are thus unfavorable 
and all nucleated bubbles with radii less than critical radius collapse and die 
out. But those with radii greater than the critical radius expand until they 
coalesce with each other. So super cooling occurs before the new phase actually 
appears and is then followed by reheating due to release of latent heat. The 
bubble nucleation rate [ [To| at temperature T is given by 

I = I e~— (1) 
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where I is the prefactor having dimension of T 4 . The prefactor used tradi- 
tionally in early Universe studies (T(| is given by I — (^)^r 4 - Csernai and 
Kapusta |ll[] have recently computed this prefactor in a coarse-grained effec- 
tive field theory approximation to QCD and give I = (^) ^ ^aw)' 2 wnere 
rj q = 14. 4T 3 is the shear viscosity in the plasma phase, £ g is a correlation length 
of order 0.7 fm in the plasma phase and Aw is the difference in the enthalpy 
densities of the two phases. In this letter we use both these prefactors for 
comparison. The critical bubble radius R c and the critical free energy W c are 
obtained by maximizing the thermodynamic work expended to create a bub- 
ble and are given by R c = Ph ( T )Z P ( T ) and W c = AnaR 2 c (T) /3 = where 
AP = Ph{T) — P q {T) is the pressure difference in hadron and quark phase. 

For simplicity we describe the quark matter by a plasma of massless u,d 
quarks and massless gluons without interaction. The long range non-perturbative 
effects are parameterized by the bag constant B. The light particles (pho- 
tons, neutrinos, and electrons) contribute equally to the pressure in both the 
phases. The pressure in the QGP phase is given by P q {T) = ^g q ^T A — B where 
g q ~ 51.25 is the effective number of degrees of freedom. In the hadronic phase 
the pressure is given by Ph{T) = ^g^j^T 4 where gt ~ 17.25, taking the three 
pions as massless. 

The fraction of the Universe h(t) which has been converted from QCD plasma 
phase to hadronic phase at the time t was first given by Guth and Weinberg fjjj 
and applied to cosmological first order phase transitions. Csernai and Kapusta 
[[n] gave a kinetic equation for calculating h(t) which takes bubble growth into 
account. If the early Universe cools to T c at time t c , then at same later time t 
the fraction of the Universe in hadronic phase is given by the kinetic equation 

h(t)= f I(T(t'))[l-h(t')}V(t',t) 

where V(t',t) is the volume of a bubble at time t which was nucleated at an 
earlier time t' and R(t) is the scale factor. This takes bubble growth into account 



R(t') 

R(t) 
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and can be given simply as 

r * 1 3 

(3) 



where v(T) is the speed of the growing bubble wall and can be taken to be 
v(T) = v a [l — 7jr]i where v Q = 3c. This has the correct behavior in that closer 
T is to T c slower do the bubbles grow. When T = |T C we have v(T) = 
the speed of sound of a massless gas. For T < |T C which occurs when there 
is large super cooling, we use the value v(T) = -4=. In this analysis collision 
and fusion of bubbles have not been taken into consideration. This seems to be 
justified as far as fusion of bubbles is concer ned. Witten JjJ and Kurki buomo 
[[To] have shown that for small enough bubbles, surface tension will cause them 
to coalesce into larger bubbles. This distance scale has been estimated to be 
given by l c — 3(^-)^( 200 ^ e y )~ 5 mm. If nucleation scale l n is larger than l c 
then we are justified in neglecting fusion of bubbles. As we will see, this is the 
case in all interesting scenarios. 

The other equation we need is the dynamical equation which couples time 
evolution of temperature to the hadron fraction h(t). We use the two Einstein's 
equations as applied to the early Universe neglecting curvature. 



R 8nG i 

i? = v— p " (4) 

R = _J_dp 

R 3w dt w 
where w = p + P is the enthalpy density of the Universe at time t. The energy 
density in the mixed phase is given by p(T) = h(t)ph(T) + [1 — h(t)]p q (T), where 
Ph and p q are the energy densities in the two phases at temperature T and 
similarly for enthalpy. We have numerically integrated the coupled dynamical 
equations (2), (4) and(5) to study the evolution of the phase transition starting 
above T c at some temperature T corresponding to time t. This time t has been 
obtained by integrating the Einstein's equation (4) and (5) and is given by 



f T dT 

t = -\l7^ZF, I , n ^ 6 ( 6 ) 
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The number density of nucleated sites at time t, is given by 

N(t) = £l(t')[l-h(t')](^ S j dt' (7) 

Therefore the typical separation between nucleation sites is l n = A^t) -1 / 3 . This 
distance scale will eventually determine the number of quarks in a nugget. This 
scale can be up to 10 12 Km depending on the parameters B and a which corre- 
spond to a distance of ~ lAMpc today. Thus the structures generated in this 
transition can be of the order of distance between galaxies today. The observable 
separation of galaxies in the Universe can be a remnant of this transition with 
the centers of the galaxies being the quark nuggets. Of course the collision of 
bubbles and their random nucleation and interaction will also lead to clustering 
of the nuggets, which can qualitatively explain the clustering of galaxies. 

To get an idea about the super cooling before nucleation begins, we can plot 
the nucleation time as a function of temperature, defined by T~^ cleation = 4 ^ c I. 
This characteristic time scale neglects bubble growth. The quark number den- 
sity in the early Universe is given by n q — -^((3)(^-)T 3 where ^ is the quark 
to photon ratio estimated from the abundance of luminous matter in the Uni- 
verse to be roughly equal to 3 x 1CU 10 . The quark nuggets which have trapped 
the quarks in the plasma phase have N q quarks at time t given by the number of 
quarks in a volume jjj^j, i-e. N q = -^y- The nucleation sites are actually ran- 
domly distributed, but we expect a distribution of quark numbers around N q . 
The average temperature at which nuggets are formed when bubbles coalesce 
is obtained by finding the average time at which the expanding bubble surfaces 
meet. Assuming a cubic lattice , this is given by setting twice the radius of the 
expanding bubbles equal to the length scale l n of the lattice spacing. We have 
done this numerically to get the corresponding time t f and temperature Tf and 
these quantities have been estimated for different values of B and a. When the 
fraction of the space occupied by the bubbles is around 50 percent, we expect 
the bubbles to meet in an ideal picture, i.e. if all bubbles are essentially nucle- 
ated at one instant which is the maximum nucleation time and they all have the 
same radius. However we have a distribution of expanding bubble sizes because 
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of the different points of time at which they were nucleated. Therefore the esti- 
mate of the time of nugget formation by treating all bubbles to be of the same 
size is an underestimate. We find that hadron fraction h(t) is only around .12 
when bubbles meet by this criteria. However we do not expect this to change 
qualitatively the broad picture of the transition and the nugget formation apart 
from reducing the formation time. 

In Fig.l we have plotted the log w of nuclcation time r as a function of 
temperature for different values of the bag pressure B, surface tension a and 
the prefactor I D . These curves clearly show when nucleation becomes large and 
how much super cooling of the Universe occurs. Fig. 2 shows the temperature 
as a function of time. It is clear from this diagram that reheating takes place 
as nucleation begins . As a increases, the super cooling is larger and the tran- 
sition takes much longer to complete with more chance of nugget formation. 
At the bottom of the curves, nucleation starts with release of latent heat and 
consequent rise in temperature. When the transition completes, the Universe 
again starts cooling. However, larger the super cooling, slower is the reheating 
as the Universe keeps on expanding. This allows the nuggets to be formed at 
a much lower temperature when bubble walls meet. From Fig. 3 we see that 
the average bubble density N(t) is initially zero and then increases with time. 
As soon as reheating starts, bubble nucleation shuts off at a particular point. 
The transition now continues only by expansion of the nucleated bubbles. The 
fall in N(t) beyond this point is due to the expansion of the Universe. Fig. 4 
shows the fraction h(t) of the Universe in hadron phase as a function of time. 
For small values of a the transition completes quickly as h(t) goes to 1. But 
for larger a it takes a larger time for h(t) to become 1. It may be compara- 
ble and even larger than the expansion time scale of the Universe. Thus the 
dilution factor is very important here as it makes the reheating much slower 
and consequently delays the transition. We also notice that in the large super 
cooling scenario the Kapusta prefactor becomes much bigger than the standard 
one by many orders of magnitude. This makes the nucleation rate as well as 
the reheating faster. In the case of low super cooling the two prefactors give al- 
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most identical results. The number of quarks in the horizon N q u at time t is 
N qH ~ n q (^Trt 3 ) ~ (^-)^^T 3 47rt 3 and we find that for all interesting cases 
N q < N qH and this number is very sensitively dependent on the surface ten- 
sion. Physically it is possible to have N q n > N q > 10 52 for some values of the 
parameters B and a. In table I below we list some physical quantities for some 
representative values of B and a. 
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Tf 


N g 


N qH 




MeV 


MeVfm- 2 


MeV 


[IS 


MeV 






m 


235 


50 


169 


12.1 


169 


2.6 x 10 28 


7 x 10 52 


8 x 10~ 3 


100 


39.5 


71.9 


2595 


13.9 


1.2 x 10 50 


3.8 x 10 56 


8.4 x 10 5 


113 


57.1 


81.3 


5138 


12.3 


9 x 10 49 


2 x 10 57 


1.7 x 10 6 



Table I. Some relevant physical quantities for some representative values of B 
and a using the standard prefactor. 

As nuggets with N q > 10 52 are expected to survive the transition, they will 
contribute to the energy density of the universe. However, nuggets of solar 
mass range need somewhat unrealistically large surface tension, specially if the 
bag pressure B is also large. 

In conclusion, detailed dynamics of the quark-hadron transition in the early 
Universe show that the evolution of the Universe does not necessarily follow 
the small super cooling scenario and certain choices of B and a can have a 
bearing on the present state of the Universe. We have explored in detail the 
possibility of nugget formation and also estimated their average separation, time 
of formation, quark content and survivability by solving a rate equation coupled 
to the Einstein's equations in an expanding universe. Clearly, the analysis 
can be improved by taking interactions into account in both the phases and 
also bubble interactions may be incorporated in the calculations. This will be 
reported elsewhere, however we believe that qualitatively the results given here 
will hold. If the nuggets studied above arc indeed formed in a much cooler 
environment, they could contribute significantly to the missing mass in the 
Universe and be candidates for dark matter. 
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Figure 1: Log of the nucleation time r in units of fm/c as a function of tempera- 
ture. Solid and dashed curves are for B 1 / 4 = WOMeV and a = 39.5MeVf m - 2 
with the standard and the Kapusta prefactors respectively. Long dashed and 
dotted curves are for B 1 ' 4 = U3MeV and a = hl.lMeVfm- 2 with the stan- 
dard and the Kapusta prefactors respectively. 
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Figure 2: Temperature as a function of time. Solid, dotted, dashed and long 
dashed curves are as in fig. 1. 
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Figure 3: Average bubble density N(t) as a function of time in units of 
(10 5 A;m 3 ) -1 . The bubble density for the standard prefactor (solid and long 
dashed curves) is normalized by multiplying with a factor of 10 3 . Curves as in 
fig. 1. 
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Figure 4: The hadron fraction as a function of time. Curves as in fig. 1. 
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